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Abstract 

The aim of this paper is first the detection of multiple abrupt changes of the 
long-range dependence (respectively self- similarity, local fractality) parameters from 
a sample of a Gaussian stationary times series (respectively time series, continuous- 
time process having stationary increments). The estimator of the m change instants 
(the number m is supposed to be known) is proved to satisfied a limit theorem 
with an explicit convergence rate. Moreover, a central limit theorem is established 
for an estimator of each long-range dependence (respectively self-similarity, local 
fractality) parameter. Finally, a goodness-of-fit test is also built in each time domain 
without change and proved to asymptotically follow a Khi-square distribution. Such 
statistics are applied to heart rate data of marathon's runners and lead to interesting 
conclusions. 

Keywords: Long-range dependent processes; Self-similar processes; Detection of abrupt 
changes; Hurst parameter; Self-similarity parameter; Wavelet analysis; Goodness-of-fit 
test. 



1 Introduction 

The content of this paper v^^as motivated by a general study of physiological signals of run- 
ners recorded during endurance races as marathons. More precisely, after different signal 
procedures for "cleaning" data, one considers the time series resulting of the evolution of 
heart rate (HR data in the sequel) during the race. The follov^^ing figure provides several 
examples of such data (recorded during Marathon of Paris 2004 by Professor V. Billat and 
her laboratory LEPHE, see jhttp:/ / www. billat . net[ ). For each runner, the periods (in ms) 
between the successive pulsations (see Fig. are recorded. The HR signal in number of 
beats per minute (bpm) is then deduced (the HR average for the whole sample is of 162 
bpm). 
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Figure 1: Heat rate signals of Athlete 1 in ms, Hertz and BPM (up), of Athletes 2, 3 and 
4 in BPM (down) 

Numerous authors have studied heartbeat time series (see for instance [24], [25j or [3j). 
A model proposed to fit these data is a trended long memory process with an estimated 
Hurst parameter close to 1 (and sometimes more than 1). In [17] three improvements 
have been proposed to such a model: 1/ data are stepped in three different stages which 
are detected using a change point's detection method (see for instance [E] or [21]). The 
main idea of the detection's method is to consider that the signal distribution depends 
on a vector of unknown characteristic parameters constituted by the mean and the vari- 
ance. The different stages (beginning, middle and end of the race) and therefore the 
different vectors of parameters, which change at two unknown instants, are estimated. 2/ 
during each stage, a time-continuous Gaussian process is proposed for modelling the de- 
trended time series. This process is a generalization of a fractional Gaussian noise (FGN) 
also called locally fractional Gaussian noise such that, roughly speaking, there exists a 
local-fractality parameter if G R (corresponding to Hurst parameter for FGN) only for 
frequencies |^| G [/™„ , fmax] with < fmin < fmax (see more details below). 3/ this 
parameter H which is very interesting for interpreting and explaining the physiological 
signal behaviours, is estimating from a wavelet analysis. Rigorous results are also proved 
providing a central limit theorem satisfied by the estimator. 

In order to improve this study of HR data and since the eventual changes of H values are 
extremely meaningful for explaining the eventual physiological changes of the athlete's 
HR during the race, the detection of abrupt change of H values is the aim of this pa- 
per. By this way the different stages detected during the race will be more relevant for 
explaining the physiological status of the athlete than stages detected from changes in 
mean or variance. For instance, the HR of a runner could decrease in mean even if the 
"fiuctuations" of the HR does not change. 

In this paper, an estimator of m instants (m G IN*) of abrupt changes of long-range 
dependence, self-similarity or local-fractality (more details about these terms will be pro- 
vided below) is developed for a sample of a Gaussian process. Roughly speaking, the 
principle of such estimator is the following: in each time's domain without change, the 
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parameter of long-range dependence (or self-similarity or local self-fractality) can be esti- 
mated from a log-log regression of wavelet coefficients' variance onto several chosen scales. 
Then a contrast defined by the sum on every m + 1 possible zones of square distances 
between points and regressions lines is minimized providing an estimator of the m instants 
of change. Under general assumptions, a limit theorem with a convergence rate satisfied 
by such an estimator is established in Theorem 12.11 

Moreover, in each estimated no-change zone, parameters of long-range dependence (or 
self-similarity or local self-similarity) can be estimated, first with an ordinary least square 
(OLS) regression, secondly with a feasible generalized least square (FGLS) regression. 
Central limit theorems are established for both these estimators (see Theorem 12.21 and 
Proposition 12.31 below) and confidence intervals can therefore be computed. The FGLS 
estimator provides two advantages: from the one hand, its asymptotic variance is smaller 
than OLS estimator one. From the other hand, it allows to construct a very simple 
(Khi-square) goodness-of-fit test based on a squared distance between points and FGLS 
regression line. The asymptotic behavior of this test is provided in Theorem 12.41 

Then, different particular cases of Gaussian processes are studied: 

1. long-range dependent processes with abrupt changes of values of LRD parameters. 
In such time series case, a semi-parametric frame is supposed (including fractional 
Gaussian noises (FGN) and Gaussian FARIMA processes) and assumptions of limit 
theorems are always satisfied with interesting convergence rates (see Corollary 13. 2p . 

2. self-similar time series with abrupt changes of values of self-similarity parameters. 
In such case, fractional Brownian motions (FBM) are only considered. Surprisingly, 
convergences of estimators are only established when the maximum of differences 
between self-similarity parameters is sufficiently small. Simulations exhibit a non 
convergence of the estimator of instant change when a difference between two pa- 
rameters is too large (see Corollary 13.411 . 

3. locally fractional Gaussian processes with abrupt changes of values of local-fractality 
parameters. In such a continuous time processes' case, a semi-parametric frame is 
supposed (including multiscale fractional Brownian motions) and assumptions of 
limit theorems are always satisfied with interesting convergence rates (see Corollary 

mil. 

The problem of change-point detection using a contrast minimization was first studied in 
the case of independent processes (see for instance Bai and Perron [5j), then for weakly 
dependent processes (see for instance Bai \4\, Lavielle [19j or Lavielle and Moulines [20j) 
and since middle of 90's in the case of processes which exhibit long-range dependance 
(see for instance Giraitis et al. [13], Kokoszka and Leipus [18] or Lavielle and Teyssiere 
|21|). Of the various approaches, some were associated with a parametric framework for 
a change points detection in mean and/or variance and others where associated with a 
non-parametric framework (typically like detecting changes in distribution or spectrum). 
To our knowledge, the semi-parametric case of abrupt change detection for long-range 
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dependent or self-similarity parameter is treated here for the first time. 
However, in the literature different authors have proposed test statistics for testing the no- 
change null hypothesis against the alternative that the long-memory parameter changes 
somewhere in the observed time series. Beran and Terrin [lO] proposed an approach 
based on the Whittle estimator, Horvath and Shao [H] obtained limit distribution of the 
test statistic based on quadratic forms and Horvath [15] suggested another test based 
on quadratic forms of Whittle estimator of long-memory parameter. The goodness-of-fit 
test presented below and which satisfies the limit theorem 12.41 also allows to test if the 
long-range memory (or self-similarity or local-fractality) parameter changes somewhere in 
the time series. 

Our approach is based on the wavelet analysis. This method applied to LRD or self- 
similar processes for respectively estimating the Hurst or self-similarity parameter was 
introduced by Flandrin [T^ and was developed by Abry, Veitch and Flandrin [2] and 
Bardet et al. [9]. The convergence of wavelet analysis estimator was studied in the case 
of a sample of FBM in [6], and in a semi-parametric frame of a general class of stationary 
Gaussian LRD processes by Moulines et al. [22] and Bardet et al. [9]. Moreover, wavelet 
based estimators are robust in case of polynomial trended processes (see Corollary 12.11) 
and is therefore very interesting for studying stochastic fiuctuations of a process without 
taking care on its smooth variations. 

A method based on wavelet analysis was also developed by Bardet and Bertrand [7] in 
the case of multiscale FBM (a generalization of the FBM for which the Hurst parameter 
depends on the frequency as a piecewise constant function) providing statistics for the 
identification (estimation and goodness-of-fit test) of such a process. Such a process was 
used for modelling biomechanics signals. In the same way, the locally fractional Gaussian 
process (a generalization of the FBM for which the Hurst parameter, called the local- 
fractality parameter, is constant in a given domain of frequencies) was studied in [T7| 
for modelling HR data during the three characteristics stages of the race. An increasing 
evolution of the local-fractality parameter during the race was generally showed for any 
runner from this method. Using the method of abrupt change detection of local-fractality 
parameter H developed in Corollary 13.61 this result is confirmed by estimations of H for 
each runner even if the change's instants seem to vary a lot depending on the fatigue of 
the runner (see the application to HR's time series in Section [3|). 

The paper is organized as follows. In Section [2], notations, assumptions and limit theo- 
rems are provided in a general frame. In Section [3l applications of the limit theorems to 
three kind of "piecewise" Gaussian process are presented with also simulations. The case 
of HR data is also treated. Section [4] is devoted to the proofs. 
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2 Main results 



2.1 Notations and assumptions 

First, a general and formal frame can be proposed. Let {Xtjt^T be a zero-mean Gaussian 
process with T = IN or T = R and assume that 

(Xq, Xs^,X2s,^, Xnsn) is known with 6n = I ov 6n — > 0, 

W— >oo 

following data are modeled with a time series (T = ISf) or a continuous time process 
T — H. In the different proposed examples X could be a stationary long memory time 
series or a self-similar or locally fractional process having stationary increments. 
For estimations using a wavelet based analysis, consider : R ^ R a function called "the 
mother wavelet". In applications, ip is a function with a compact (for instance Daubeshies 
wavelets) or an essentially compact support (for instance Lemarie- Meyer wavelets). For 
{Xt)teT and (a, b) e R^ x R, the wavelet coefficient of X for the scale a and the shift b is 

dx{a,b) / ,l^{^-^)X{t)dt. 

When only a discretized path of X is available (or when T = IN), approximations ex (a, b) 
of dx{(i,b) are only computable. We have chosen to consider for (a, 6) G M*^ x 

ex(a,6) := ^ (1) 
Va a 

which is the formula of wavelet coefficients computed from Mallat's algorithm for com- 
pactly supported discrete (a e 2"^) wavelet transform (for instance Daubeshies wavelets) 
when N is large enough and nearly this formula for discrete wavelet transform with an 
essentially compact support (for instance Lemarie-Meyer wavelets). Now assume that 
there exist m e IN (the number of abrupt changes) and 

• — Tq < T* < . . . < < r^+i = 1 (unknown parameters) ; 

• two families {c(j)o<j<m G R™"*"^ and (/3*)o<j<m G (0,00)™+^ (unknown parameters); 

• a sequence of "scales" (an)neiN £ (known sequence) satisfying a„ > amin for all 
n e IN, with amin > 0, 

such that for all j = 0, 1, . . . , m and k e -D^(j) C [[NSnt*] , [NSntJ^i]] , 

E[e^(ajv, k)] ~ p* ■ (ajv)"^* when AT ^ 00 and NSn ^ 00. (2) 

Roughly speaking, for N & ¥1* the change instants are [N6nt*] for j = l,...,m, the 
variance of wavelet coefficients follows a power law of the scale, and this power law is 
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piecewise varying following the shift. Thus piecewise sample variances can be appropriated 
estimators of parameters of these power laws. Hence let us define 

[fc'/ajvl-l 

•= Yl 4K,a^p) for < A; < A;' < iV5^. (3) 

p=[fc/ajv] 



Now set < ri < . . . < r£ with ^ G IN* and let us suppose that a multidimensional central 

'l<i<£' 



limit theorem can also be established for (^^'(rj a7v))j t-e. 



with [NSnt*] < k < k' < [N5nt;+;\ and it exists T^^\a*,n, . . . ,re) = {7pg)i<p^q<i a 
(i X i) matrix not depending on such that a t-^ T^^\a,ri, . . . ,ri) is a continuous 
function, a positive matrix for all a and 



{el''\k,k'))^^^^^ ^^X{0,T^^\a*,n,...,n)) when k' - k ^ oo. (5) 
With the usual Delta-Method, relation (Jlj) implies that for 1 < z < £, 



log (5f (r, a^)) = log(/?;) + a* log (r, a^) + ef\k, k'), (6) 

for [N6nt*] < k < k' < [N6nt*^i] and the limit theorem ([5]) also holds. This is a linear 
model and therefore a log-log regression of (S^' {ri Oiv)) . onto (rj oat) . provides an estima- 
tor of a* and log(/?*). 

The first aim of this paper is the estimation of unknown parameters (T*)j, and 
Therefore, define a contrast function 

m £ 

f^Jv((«j)o<i<m, iPj)o<j<m, {kj)l<j<m) = ( log (5"^^'^' (r^ Oiv)) - log(ri Oat) +log 

[ • (a,)o<,<,n G C R'^+i 

with <^ • (/?,)o<,<m G 5"'+' C (0, oo)'"+i 

The vector of estimated parameters aj, (3j and kj (and therefore Tj) is the vector which 
minimizes this contrast function, «. e., 

(('^i)o<j<m, (/5j)o<j<m, (^j)l<j<m) 

:= Argmin|?77v((aj)o<j<m, (/5i)o<j<m, (fci)i<i<m)} in x 5"^+^ x ir„(Arj(7) 

Tj- := kj/{N5N) for 1 < J < m. (8) 
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For a given {kj)i<j<m, it is obvious that (5?j)o<i<m and (log/5j)o<j<m are obtained from a 
log-log regression of {Sf.^.'*'^ {rittN))- onto (rjOAr)., i.e. 



. Therefore the 



/ log(riaAr) 1 ^ 
with Y^^+' := (log {S^'.^^ri ■ aN)))^^-^^ and L^^ := ; : 

y log(r^a7v) 1 J 

estimator of the vector (fcj)i<j<m is obtained from the minimization of the contrast 

GN^ki, k2, . . . , km) '■= UN{{o:j)o<j<m, {Pj)o<j<m, {kj)i<j<m) (9) 

=^ (kj)i<j<m = Argmin|GAr(A;i, k2, ■ ■ ■ , k^), {kj)i<j<m e Km{N)^. (10) 

2.2 Estimation of abrupt change time- instants (T*)i<j<m 

In this paper, parameters (a*) are supposed to satisfied abrupt changes. Such an hypoth- 
esis is provided by the following assumption: 

Assumption C: Parameters (a*) are such that la*^,^ — a;*| 7^ for all j = 0, 1, . . . , m — 1. 
Now let us define: 

r* := (ti, . . . ,r^), r := (n, . . . , f^) and ||r||„ := max (|ri|, . . . , Ix^l). 
Then r converges in probability to r* and more precisely, 

Theorem 2.1 Let £ G IN \ {0,1,2}. // Assumption C and relations ^ and ^ 
hold with (a*)o<j<m such that a* G [a, a'] and a < a' for all j = 0,...,m, then if 
^i+2ia -a)^j^ ^^yi — ^ 0^ for all {vn)n satisfying Wat ■ a^"^^" ~"^(A^ ^at)"^ — > 0, 

P(vN\\r* -I.\\m>ri) — ^0 forallr]>0. (11) 
Several examples of applications of this theorem will be seen in Section [3l 

2.3 Estimation of parameters {aj)o<j<m and {Pj)o<j<m 

For j = 0,1,..., m, the log-log regression of (5'|-'+'('^i«A'))i<i<£ onto {riaN)i<i<e pro- 
vides the estimators of a* and (3*. However, even if tj converges to r*, kj = N5n ■ tj 
does not converge to k* (except if iV = o{vn) which is quite impossible), and therefore 



7 



P([/cj, kjj^i] C [k*, A^^+i]) does not tend to 1. So, for j = 0, 1, . . . , m, define kj and k'j such 
that 

L=L + ^ and ~E = kj+^-^ =^ F([L,k]c[klk*]) ^1, 

from m with 7^ = 1/2. Let 6* := j^^^^, ^ and 6, := (U,- L,y'lJ,-Y^^ :^ ^ 
Thus, estimators ('5j)o<j<m and (/3j)o<j<m satisfy 



Theorem 2.2 Under the same assumptions as in Theorem \2.1\. for j = 0, . . . ,m 



■(®^-®^*) ^-^Ar(0,S(^)(a*,ri,...,r,)) (12) 



Oat 

With S(^)(a*,ri, . . . ,r,) := (L; ■ U)-^L\ ■ T^^\a], n, . . . ,r,) ■ Li ■ {L[ ■ L^)-\ 

A second estimator of 0* can be obtained from feasible generalized least squares (FGLS) . 
Indeed, the asymptotic covariance matrix T^^\a*, ri, . . . , r^) can be estimated with the ma- 
trix f*^-'^ := T^^\aj, ri, . . . , r^) and f*^-'^ T''^\a*^ ri, . . . , r^) since a ^ r'^-'^(a, ri, . . . , r^) 

N~*oo ■' 

is supposed to be a continuous function and otj a*. Since also a ^— T^^^a, ri, . . . , r^) 

is supposed to be a positive matrix for all a then 

Then, the FGLS estimator Qj of G* is defined from the minimization for all B of the 



following criterion 

II 1^ - L., . e 11^,,= (i| - L., . e)' . (f (^))-^ . (1^^ - L., . e). 

and therefore 



Proposition 2.3 Under the same assumptions as in Theorem \2.Sl for j = 0, . . . ,m 

„-!i^(M/»K'.n... ...,)) (13) 

with M^^\a*, n, . . . , r,) := {L[ ■ {T^^\a*, ri, . . . , r,))"' • Li)"' < n, . . . , r,) (with 

order's relation between positive symmetric matrix). 

Therefore, the estimator Qj converges asymptotically faster than G^; cxj is more interesting 
than ttj for estimating a* when N is large enough. Moreover, confidence intervals can be 
easily deduced for both the estimators of 0*. 
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2.4 Goodness-of-fit test 

For j = 0, . . . , m, let T^^^ be the FGLS distance between both the estimators of i.e. 

the FGLS distance between points ( log(rj otv), log (5'-'') ) and the FGLS regression 

V '^j / i<i<e 

line. The following limit theorem can be established: 

Theorem 2.4 Under the same assumptions as in Theorem \2.1l for j = 0, . . . ,m 

T« = '""^'t''"^' II ^'^ - ^» ■ ^' - 2). (14) 

Mutatis mutandis, proofs of Proposition 12.31 and Theorem 12.41 are the same as the proof 
of Proposition 5 in |7j. This test can be applied to each segment [kj, kj[. However, under 
the assumptions, it is not possible to prove that a test based on the sum of T^^^ for 
j = 0, . . . ,m converges to a + 1)(^ ~ 2)) distribution (indeed, nothing is known 

about the eventual correlation of (Y-^)^,., ). 

V kj ' 0<j<m-' 

2.5 Cases of polynomial trended processes 

Wavelet based estimators are also known to be robust to smooth trends (see for instance 
[T]). More precisely, assume now that one considers the process Y = {Yt, t E T} satisfying 
Yf = Xt + P{t) for all t e T where P is an unknown polynomial function of degree p G IN. 
Then, 

Corollary 2.1 Under the same assumptions as in Theorem \2. 1\ for the process X , and if 
the mother wavelet ip is such that J f il){t)dt = for r = 0, 1, . . . then limit theorems 
and ^ hold for X and for Y. 

Let us remark that Lemarie-Meyer wavelet is such that / 1^ ip(t)dt = for all r G IN. 
Therefore, even if the degree p is unknown, Corollary 12.11 can be applied. It is such the 
case for locally fractional Brownian motions and applications to heartbeat time series. 



3 Applications 

In this section, applications of the limit theorems to three kinds of piecewise Gaussian 
processes and HR data are studied. Several simulations for each kind of process are 
presented. In each case estimators (Tj)j and {aj)j are computed. To avoid an overload 
of results, FGLS estimators which are proved to be a little more accurate than 

{(yj)j are only presented in one case (see Table [2]) because the results for (aj)j are very 
similar to {<yj)j ones but are much more time consuming. For the choice of the number 
of scales £, we have chosen a number proportional to the length of data (0.15 percent of 
which seems to be optimal from numerical simulations) except in two cases (the case 
of goodness-of-fit test simulations for piecewise fractional Gaussian noise and the case 
of HR data, for which the length of data and the employed wavelet are too much time 
consuming). 
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3.1 Detection of change for Gaussian piecewise long memory pro- 
cesses 

In the sequel the process X is supposed to be a piecewise long range dependence time 
series (and therefore 6^ = 1 for all N G IN). First, some notations have to be provided. 
For Y = {Yt)teK a Gaussian zero mean stationary process, with r(t) = Fi(Yo-Yt) for t G IN, 
denote (when it exists) the spectral density / of F by 

/(A) = 7^- Vr(fc)-e-^'=^ for AG Ac [-TT.vr]. 

In the sequel, the spectral density of Y is supposed to satisfy the asymptotic property, 

/(A) ~ C ■ when A 0, 
A 

with C > and D G (0, 1). Then the process Y is said to be a long memory process 
and its Hurst parameter is if = (1 + D)/2. More precisely the following semi-parametric 
framework will be considered: 

Assumption LRD(D): F is a zero mean stationary Gaussian process with spectral 
density satisfying 

/(A) = |A|-^-/*(A) for all A G [-tt, 0[U]0, tt], 

with /*(0) > and /* is such that |/*(A) - /*(0)| < C2 ■ |A|2 for all A G [-tt, vr] with 
C2 > 0. 

Such assumption has been considered in numerous previous works concerning the esti- 
mation of the long range parameter in a semi-parametric framework (see for instance 
Robinson, 1995,, Giraitis et al, 1997, Moulines and Soulier, 2003). First and famous ex- 
amples of processes satisfying Assumption LRD(D) are fractional Gaussian noises (FGN) 
constituted by the increments of the fractional Brownian motion process (FBM) and the 
fractionally autoregressive integrated moving average FARIMA[p, c?, g] (see more details 
and examples in Doukhan et al. [TT]). 

In this section, X = {Xt)tem is supposed to be a Gaussian piecewise long-range dependent 
process, i.e. 

• there exists a family (-D*)o<j<m G (0, 1)™+^; 

. for all J = 0, . . . ,m, for all k G {[iVr/], [Nr*] + 1, . . . , [Nr*^,] - l}, X^ = X^i\^^., 
and = (X^^^)tg]N satisfies Assumption LRD(L'*). 

Several authors have studied the semi-parametric estimation of the parameter D using 
a wavelet analysis. This method has been numerically developed by Abry et al. (1998, 
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2003) and Veitch et al. (2004) and asymptotic results are provided in Bardet et al. (2000) 
and recently in Moulines et al. (2007) and Bardet et al. (2007). The following results 
have been developed in this last paper. The "mother" wavelet ip is supposed to satisfy 
the following assumption: first ip is included in a Sobolev space and secondly ip satisfies 
the admissibility condition. 

Assumption Wi : ^p : E t-^ R with [0, l]-support with ^p{0) = ^p{l) = and ip{t) dt = 
and such that there exists a sequence {ipe)£(zx satisfying ip{X) = J2eez ^ ]L^([0, 1]) 

and E.ez(l + l^lf/'l^^l <oo- 

For ease of writing, ip is supposed to be supported in [0, 1]. By an easy extension the 
following propositions are still true for any compactly supported wavelets. For instance, 
ip can be a dilated Daubechies "mother" wavelet of order d with d > 6 to ensure the 
smoothness of the function ip. However, the following proposition could also be extended 
for "essentially" compactly supported "mother" wavelet like Lemarie-Meyer wavelet. Re- 
mark that it is not necessary to choose ip being a "mother" wavelet associated to a 
multi-resolution analysis of IL^(IR) like in the recent paper of Moulines et al. (2007). The 
whole theory can be developed without resorting to this assumption. The choice of ip is 
then very large. Then, in Bardet et al. (2007), it was established: 

Proposition 3.1 Let X be a Gaussian piecewise long-range dependent process defined as 
above and (an)ne]N be such that N/a^ — > oo and ■ N~^/^ — > oo. Under Assump- 



tion Wi, limit theorems ^ and ^ hold with a* = D* and (3* = log ( /*(0) |^/'(n) | 



\u\ 



^duj for all j = 0,1, ... ,m and with dpg = GCD{r,p , r^) for all {p, g) G {1, . . . 
2(rprg)^"-°j* >^ / ip{urp)ilj{urq) u^^j cos{udpqm)du 



'PI fj 

Ujpq 



m=—oo 



As a consequence, the results of Section [2] can be applied to Gaussian piecewise long-range 
dependent processes: 

Corollary 3.2 Under assumptions of Proposition \3. il and Assumption C, for all < k. < 
2/15, ifaN = N^+^/^ and Vn = N^'^-^^ then (EP, HE), HE) and g hold. 

Thus, the rate of convergence of r to r* (in probability) is A^2/5-3k for < k as small as 
one wants. Estimators Dj and Dj converge to the parameters D* following a central limit 
theorem with a rate of convergence A^/^"''/^ for < as small as one wants. 

Results of simulations: The following Table [T] represents the change point and parameter 
estimations in the case of a piecewise FGN with one abrupt change point. We observe 
the good consistence property of the estimators. Kolmogorov-Smirnov tests applied to 
the sample of estimated parameters lead to the following results: 
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1. the estimator ri can not be modeled with a Gaussian distribution; 

2. the estimator Hj seems to follow a Gaussian distribution. 



N = 20000, Ti = 0.75, Do = 0.2 and = 0.8 


fi ar, VMSE 


Do do, VMSE 


Di do, VMSE 


0.7605 0.0437 0.0450 


0.2131 0.0513 0.0529 


0.7884 0.0866 0.0874 



Table 1: Estimation of ri. Do and Di in the case of a piecewise FGN {Ho = 0.6 and 
Hi = 0.9) with one change point when = 20000 and £ = 30 (50 realizations) 

The distribution of the test statistics T(°) and T^^) (in this case £ = 20 and = 20000 
and 50 realizations) are compared with a Chi-squared-distribution with eighteen degrees of 
freedom. The goodness-of-fit Kolmogorov-Smirnov test for T^^^ to the x^(18)-distribution 
is accepted (with 0.3459 for the sample of T^^^ and p = 0.2461 for T^^)). In this case and 
for the same parameters as in Table [U the estimator Dj seems to be a little more accurate 
than Dj (see Table ED . 



fi dr, VMSE 


Do aDo VMSE 


Di an, VMSE 


0.7652 0.0492 0.0515 


0.1815 0.0452 0.0488 


0.8019 0.0721 0.0722 



Table 2: Estimation of Do and Di in the case of a piecewise FGN [Do = 0.2 and Di = 0.8) 
with one change point when N = 20000 and i = 20 (50 realizations) 

Simulations are also applied to a piecewise simulated FARIMA(0,(ij,0) processes and re- 
sults are similar (see Table [3]). The following Figure [2] represents the change point instant 
and its estimation for such a process with one abrupt change point. 




Figure 2: Detection of the change point in piecewise FARIMA(0,(ij,0) (for the first segment 
do = 0.1 (Do = 0.2) for the second di = 0.4 (D^ = 0.8)) 
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N = 20000, Ti = 0.75, Do = 0.2 and Di = 0.8 


fi dr, VMSE 


Do do, VMSE 


Di do, VMSE 


0.7540 0.0215 0.0218 


0.1902 0.0489 0.0499 


0.7926 0.0761 0.0764 



Table 3: Estimation of ri. Do and Di in the case of piecewise FARIMA(0,(ij,0) {do = 0.1 
and di = 0.4) with one change point when = 20000 and £ — 30 (50 reahzations) 

3.2 Detection of abrupt change for piecewise Gaussian self-similar 
processes 

Let us recall that = {B^)t£u is a fractional Brownian motion (FBM) with two param- 
eters H e (0, 1) and > when Bh is a Gaussian process having stationary increments 
and such as 

Var(Sf)) =(72|i|2^ \/t e E. 

It can be proved that Bh is the only Gaussian self-similar process having stationary 
increments and its self-similar parameter is H (a process Y = (Yt)t^E is said to be a 
i^s-self-similar process if for all c > and for all (ti, . . . ,tk) G E^ where A; G IN*, the 
vector (l^ti, . . . , yctk) has the same distribution than the vector (F^^, . . . , 1^^.)). 

Now, X will be called a piecewise fractional Brownian motion if: 

• there exist two families of parameters {H*)o<j<m G (0, 1)™+^ and {cr'f')o<j<m G 
(0,oo)"^+^ 

. for all J = 0, ... , m, for all t G [[Nr*], [Nr*] + 1, . . . , [Nr*^^] - l] , X, = x'^^^^^ 
and X^-^) = {X[^^)te-R is a FBM with parameters H* and af. 

The wavelet analysis of FBM has been first studied by Flandrin (1992) and developed 
by Abry (1998) and Bardet (2002). Following this last paper, the mother wavelet -0 is 
supposed to satisfy: 

Assumption W'2,: : R — > R is a piecewise continuous and left (or right)-differentiable 
in [0,1], such that \%l)'{t~)\ is Riemann integrable in [0,1] with ip'{t~) the left-derivative 
of -0 in t, with support included in [0, 1] and J^tPilj{t) dt — t''''il){t) dt — for p — 0, 1. 

As in Assumption Wi, ip is supposed to be supported in [0, 1] but the following proposi- 
tions are still true for any compactly supported wavelets. Assumption W2 is clearly weaker 
than Assumption Wi concerning the regularity of the mother wavelet. For instance, ijj 
can be a Daubechies wavelet of order d with li > 3 (the Haar wavelet, i.e. d = 2, does not 
satisfy j^tip{t)dt = 0). Another choice could be infinite support wavelets with compact 
effective support (it is such the case with Meyer or Mexican Hat wavelets) but the proof 
of the following property has to be completed. 
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Proposition 3.3 Assume that X is a piecewise FBM as it is defined above and let 
(Xi, X2, . . . , Xtv) he a sample of a path of X (therefore 6n = Ij- Under Assumption 
W2, if (an)new is such that N/a^ — > c)o and ■ N~'^^^ — > 00, then limit theorems 

N-*oo N^oo 



(gP, ^ and (M) hold with a* = 2H* + 1 and P* = log {-^ ij{t)ij{t')\t-t'\^^Ut dt') 
for all j = 0,1, ... ,m and with dpq = GCDijp , r^) for all {p,q) G {1, . . . , d}, 

(j) _ 2(ipg / j'l j'l i){t)%l){t') \kdpg + Tpt - rqt'f"^ dt dt' 



fc=— oo 



r^^{t)^{t')\t-tf''Utdt' 



Then, Theorem 12.11 can be applied to piecewise FBM but 2(a' — a) + 1 = 2(supj a 
inf, a*) + 1 has to be smaller than 3 since ■ N^^^^ — > oo. Thus, 



Corollary 3.4 Let A := | sup^ H* — infj H* \ . If A < 1/2, under assumptions of Propo- 
sition \3.3\ and Assumption C, for all < k < j+aa ~ 1> '^n = N'^^'^'^'^ and vn = 

p^2/3{l~2A)-.(2+4A) ^^^^ ^J^Jj, ^^ , ^^ and g hold. 

Thus, the rate of convergence of r to r* (in probability) can be iv2/3(i-2yi)-ft' < 
as small as one wants when = A^V3+k'/(2+4A)_ 

Remark: This result of Corollary 13.41 is quite surprising: the smaller A, i.e. the smaller 
the differences between the parameters Hj, the faster the convergence rates of estimators 
Tj to T*. And if the difference between two successive parameters Hj is too large, the 
estimators Tj do not seem to converge. Following simulations in Table [5] will exhibit this 
paroxysm. This induces a limitation of the estimators' using especially for applying them 
to real data (for which a priori knowledge is not available about the values of H*). 

Estimators Hj and Hj converge to the parameters H* following a central limit theo- 
rem with a rate of convergence N^/^~'^/'^ for < k as small as one wants. 

Results of simulations: The following Table [4] represent the change point and parame- 
ter estimations in the case of piecewise FBM with one abrupt change point. Estimators 
of the change points and parameters seem to converge since their mean square errors 
clearly decrease when we double the number of observations. 

For testing if the estimated parameters follow a Gaussian distribution, Kolmogorov- 
Smirnov goodness-of-fit tests (in the case with N = 10000 and 50 replications) are applied: 

1. this test for Hq is accepted as well as for Hi and the following Figure [3] represents 
the relating distribution. 

2. this is not such the case for the change point estimator Ti for which the hypothesis 
of a possible fit with a Gaussian distribution is rejected {KStest = 0.2409) as showed 
in the Figure [3] below which represents the empirical distribution function with the 
correspondant Gaussian cumulative distribution function. 
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N = 5000 


N = 10000 




n 




VMSE 


n 




VMSE 


0.4 


0.4467 


0.0701 


0.0843 


0.4368 


0.0319 


0.0487 




Ho 




y/MSE 


Ho 




VMSE 


0.4 


0.3147 


0.0404 


0.0943 


0.3761 


0.0452 


0.0511 


Hi 


Hi 




VMSE 


Hi 




VMSE 


0.8 


0.7637 


0.0534 


0.0645 


0.7928 


0.0329 


0.0337 



Table 4: Estimation of ri, Ho and Hi in the case of piecewise FBM with one change point 
when = 5000 (100 reahzations) and = 10000 (50 reahzations) 

12r Empirical CDF 

_ 0.9 - /y 









m 

\ 


\ 






1 













Figure 3: Left: Modehng of sample estimations of Ho with normal distribution; Right: 
Comparison of the generated empirical cumulative distribution for ti (when N= 10000) 
and the theoretical normal distribution. 



From the following example in Table [5l we remark that the estimated parameters seem 
to be non convergent when the difference between the parameters Hj is too large. 



N = 5000, Ti = 0.6, Ho = 0.1 and Hi = 0.9 


fi dr, VMSE 


Ho dn, VMSE 


Hi an, VMSE 


0.5950 0.1866 0.1866 


-0.1335 0.0226 0.2346 


0.6268 0.4061 0.4894 



Table 5: Estimation of ri. Ho and Hi (when Hi —Ho = 0.8 > 1/2) in the case of piecewise 
FBM with one change point when = 5000 (50 realizations) 

Simulations for goodness-of-fit tests T^^^ provide the following results: when = 5000, 
the drawn distributions of the computed test statistics (see Figured!) exhibit a Khi-square 
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distributed values (x^(5) since e = 7) and 95% of the 100 of the values of T^o) and T^^) do 
not exceed X95%(5) = 11.0705. These results are also validated with Kolmogorov-Smirnov 
tests. 




Figure 4: Testing for distribution in the first detected zone (left) and the second 

detected zone (right) (50 realizations when = 5000) 

The results below in Table [6] are obtained with piecewise fractional Brownian motion 
when two change points are considered. As previously, both the KStest tests for deciding 
whether or not samples of both estimated change points is consistent with Gaussian 
distributions are rejected. However, such KStest tests are accepted for Hj samples. A 
graphical representation of the change point detection method applied to a piecewise 
FBM is given in Figure O 





= 5000 


N = 10000 








VMSE 






VMSE 


0.3 


0.3465 


0.1212 


0.1298 


0.3086 


0.0893 


0.0897 








VMSE 


% 




VMSE 


0.78 


0.7942 


0.1322 


0.1330 


0.7669 


0.0675 


0.0687 




Ho 




VMSE 






VMSE 


0.6 


0.5578 


0.0595 


0.0730 


0.5597 


0.0449 


0.0604 


Hi 


Hi 




VMSE 


Hi 




VMSE 


0.8 


0.7272 


0.0837 


0.1110 


0.7633 


0.0813 


0.0892 




H2 




VMSE 


H2 




VMSE 


0.5 


0.4395 


0.0643 


0.0883 


0.4993 


0.0780 


0.0780 



Table 6: Estimation of ri, T2, Hq, Hi and H2 in the case of piecewise FBM with two 
change points when N = 5000 and = 10000 (50 realizations) 

The distribution of the test statistics T(°\ T^^) and T^^) (in this case ^ = 10, = 10000 
and 50 realizations) are compared with a Chi-squared-distribution with eight degrees of 
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1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 



Figure 5: (left)Detection of the change point in piecewice FBM{Hj) (ri = 0.3, T2 = 0.78, 
Hq = 0.6, Hi = 0.8 and H2 = 0.5). The change points estimators are ri = 0.32 and T2 = 
0.77. (right) Representation of log-log regression of the variance of wavelet coefficients 
on the chosen scales for the three segments {Hq = 0.5608 (*), Hi = 0.7814 (<l) and 
H2 = 0.4751 (o)) 

freedom. The goodness-of-fit Kolmogorov-Smirnov test for T^^^ to the x^(8)-distribution 
is accepted (with p = 0.4073 for the sample of T^o), p = 0.2823 for T^^) and p = 0.0619 

for r(2)). 

3.3 Detection of abrupt change for piecewise locally fractional 
Gaussian processes 

In this section, a continuous-time process X is supposed to model data. Therefore assume 
that {Xsj^, X2Sf^, ■ ■ ■ , Sn) is known, with 6n — > and N6n — ^ 00. A piecewise 

locally fractional Gaussian process X = {Xt)t£R+ is defined by 

Xt:= [ ^^^^i^W{dO for te[T*NSN,r;^iNSN) (15) 

where the functions pj : R [0, 00) are even Borelian functions such that for all j = 
0, 1, . . . ,m,: 

• pAO = A lel^'^+'Z' for lei e [frmn, fma.] with H* G R, a* > 0; 
a. 

and W{dx) is a Brownian measure and W{d$,) its Fourier transform in the distribution 
meaning. Remark that parameters H*, called local-fractality parameters, can be sup- 
posed to be included in R instead the usual interval (0,1). Here < fmin < fmax are 
supposed to be known parameters. Roughly speaking, a locally fractional Gaussian pro- 
cess is nearly a self-similar Gaussian process for scales (or frequencies) included in a band 
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of scales (frequencies). 

For locally fractional Gaussian process already studied in Bardet and Bertrand (2007) 
and Kammoun et al. (2007), the mother wavelet is supposed to satisfy 

Assumption W^: ^Ij : R E is a C°°(R) function such that for all m G IN, /^^ |t™?/'(t) | dt < 
oo and the Fourier transform ^ of is an even function compactly supported on [— /x, — A]U 
[A, yu] with < A < /i. 

These conditions are sufficiently mild and are satisfied in particular by the Lemarie-Meyer 
"mother" wavelet. The admissibility property, i.e. f^ip{t)dt = 0, is a consequence of the 
second condition and more generally, for all m G IN, J^t"^ilj{t)dt = 0. 

Since the function ifj is not a compactly supported mother wavelet, wavelet coefficients 
dx{a,b) can not be well approximated by ex{a,b) when the shift b is close to or NSn- 
Then, a restriction (a^) of sample wavelet coefficient's variance (ajy) has to be 
defined: 

[{k'~w{k'~k))/aM]-l 

Ski(^N) ■■= n - 2w)k' - k ^ (^x{(^N, aNp) with 0<w < 1/2. 

^ p=[{k+w{k'-k))/aN]+w 



Proposition 3.5 Assume that X is a piecewise locally fractional Gaussian process as it is 
defined above and {Xsj^, X2Sj^, Xnsj^) is known, with N{6n)^ — > 0andN6N — ^ oo. 

Af^oo N^oo 

Under Assumptions W3 and C, using S^' (a^) instead of Sj^'{aN), if x < (^"i^d 
^. ^ frr^ I (^frn^ _ for i = 1, . . . J with ttj^ = 1 for all G IN, then limit theorems 

(gP, ^ and ^ hold with a* = 2H* + 1 and (3* = \og[-^J^ 1^(^)1^ \u\-^'^"^du) for 
all j = 0, 1, . . . , m, for all (p, g) G {1, . . . , i}, 

{l-2w){r,r,rUu^ JJ^iu)\'\u\-'-'^^du ) ^ ^ 

Theorem 12.11 can be applied to a piecewise locally fractional Gaussian process without 
conditions on parameters H*. Thus, 

Corollary 3.6 Under assumptions of Proposition 1 3'. 51 and Assumption C, then for all 
0<K<l,if5N = Ar-i/2-« and vn = A^^/^-k ^^^^ ^j^]), ^M), ^W and g hold. 

Therefore the convergence rate of r to r* (in probability) is as well close to A^^^^ as one 
wants. Estimators Hj and Hj converge to the parameters H* following a central limit 
theorem with a rate of convergence N^/^~'^/'^ for < k as small as one wants. 
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3.4 Application to heart rate's time series 

The study of the regularity of physiological data and in particular the heartbeat signals 
have received much attention by several authors (see for instance [24], [25] or [3]). These 
authors studied HR series for healthy subjects and subjects with heart disease. In [T7] . 
a piecewise locally fractional Brownian motion is studied for modeling the cumulative 
HR data during three typical phases (estimated from Lavielle's algorithm) of the race 
(beginning, middle and end). The local-fractality parameters are estimated with wavelet 
analysis. The conclusions obtained are relatively close to those obtained by Peng, et al. 
Indeed we remarked that the local-fractality parameter increases thought the race phases 
which may be explained with fatigue appearing during the last phase of the marathon. 
In this paper, we tray to unveil in which instants the behaviour of HR data changes. The 
following Table [7] presents the results for the detection of one change point. 







Ho 




Hi 






Athl 


0.0510 


0.7880 


1 


.2376 


1.0184 


1.0562 


Ath2 


0.4430 


1.3470 


1 


.4368 


5.0644 


1.5268 


Ath3 


0.6697 


0.9542 


1 


.2182 


0.7836 


0.9948 


Ath4 


0.4856 


1.1883 


1 


.2200 


2.8966 


1.2774 


Ath5 


0.8715 


1.1512 


1 


.3014 


0.7838 


0.8748 


Ath6 


0.5738 


1.1333 


1 


.1941 


2.2042 


0.7464 


Ath7 


0.3423 


1.1905 


1 


.1829 


0.4120 


1.5598 


Ath8 


0.8476 


1.0222 


1 


.2663 


3.1704 


0.5150 


Ath9 


0.7631 


1.4388 


1 


.3845 


9.6574 


0.5714 



Table 7: Estimated change points ri, parameters i^o, Hi and goodness-of-fit test statistics 
(T(°) for the first zone and T^^^ for the second one) in the case of one change point observed 
in HR series of different athletes. 

It is noticed that the estimator of the local-fractality parameter is generally larger on 
the second zone than on the first although the detected change point differs from an 
athlete to another (only the case of Athlete 1 seems not to be relevant). This result is 
very interesting and confirms our conclusions in [17]. Whatever is the position of change 
point, the estimation of the local-fractality parameter is larger in the second segment than 
in the first segment (see the example of HR data recorded for one athlete in Figure Ej). 

In general, the goodness-of-fit tests, with values T^^^ and T^^\ are less than ^95% (4) = 
9.4877 (except T^^^ for Ath9) when £ = 6. So, the HR data trajectory in the both zones 
seems to be correctly modeled with a stationary locally fractional Gaussian trajectory. 
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Figure 6: Evolution of local-fractality parameter estimators (observed for HR series of 
one athlete) in the two zones when the change point varies in time. 

4 Proofs 

Before establishing the proof of Theorem 12.11 an important lemma can be stated: 

Lemma 4.1 Let A; G IN \ {0, 1}, (7i)i<i<fc G (0, oo)'^ and ai > a2 > ■ ■ ■ > ak be k ordered 
real numbers. For G R^, consider the function fa^p : x G R i— R such that 

k 

fa,f}{x) := a X + /3 - log ^ 7, exp (a^ x) j for x G R. 

q=l 

Let < ti < ■ ■ ■ < tf with £ G IN \ {0, 1, 2} and (M„)„eM be a sequence of real numbers 
such that there exists m G R satisfying Un> m for all n G IN. Then there exists C > 
not depending on n such that 

inf V |/„,^(logK) > C min (l, 

1=1 



Proof of Lemma For all G R , the function fa^p is a C°°(R) function and 

, , Eg;iSg7q+i(«5-Q;g+i)^exp((ag + ag+i)x) 

= 77:^ — < °- 

(^E9=i7gexp (agx) j 

Therefore the function fa^p is a concave function such that sup(Q,^)gj^2 ^/q,/3(x) < (not 
depending on a and /?) and for all (a,/3) G R^, fa,p vanishes in 2 points at most. Thus, 
since £ > 3 and (x + tj) . are distinct points, for all x G R, it exists C(x) > not depending 
on a and (3 such that 

I 

1=1 

Therefore, since for all M > 0, 

inf I inf y"\faa{x + t,)?\ > inf {C(x)| > 0. (17) 
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Moreover, if Un — > +00, 

k 



log ( 7g exp (og log(M„))^ = log (ji exp [ai log(M„)) + 72 exp (02 log(M„)) (l + o{l)] 

q=l 

= log(7i) + ai log(M„) + 72 exp ((^2 - Oi) log(M„)) (l + o(l)). 

Thus, for n large enough, 

k 

log($^7,exp(«, logK))) -log(7i) + «ilogK)| <2j2K'"'\ (18) 



1 

2 72«r"^< 



Therefore, for all G 

|2 



fa,f3{\0g{Un)+ti)\ = \ fa^,logiyi) {^Og{Un) +U)\ + (log(7i) - /5) + («!-«) (log(M„) + tj) 
- 2/ai,log(7i)(log("n) ^ ((l0g(7l) - P) + ("1 -a)(log(M„) + ti) ) . 



Using inequalities (UHl), ^7' < |/a„iog(70 ( logK) + 1.) |' < 472<(°^-"i) and for 
all (a,/?) G R^ lim /ai,iog{7i)(log(M„) x ( (log(7i) -/?) + (ai - a)(log(M„) ) = 0. 

Then, for all (a,/?) 7^ (ai, log(7i)), lim |/a^/3(log('U„) + tj) I = 00. Consequently, for n 
large enough. 



. ,X]|.^°'/3(logK)+ti)|^ > - V |/ai,log(7i)(logK) 



which combined with (ITTll achieves the proof. □ 



Proof of Theorem \2.1[ Let = — — , k* = \N6nt*] for j = 1, . . . , m and 

vn ^ ^ 

^vy^N = {ikj)i<j<m, max \kj - k*\ > r]WN}- 

jGl,...,m ■' 



Then, for N6n large enough 

■N6_ 



pf — \\z* ~ z\\m > v) — P( max \kj — k*\ > tiwn) 

V Wn ' i6l,.--,m 

= P( min G'jv((^i)i<j<m) < min G'Ar((A;j)i<j<^) j 
< min G^((A;,)i<,<^) <G^((A:;)i<,<^)). (19) 

For j = {0, . . . , m} and = /cq < /ci < • • • < /^m < /cm+i = ^^n-, let 
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kj 



{\og{sl]^\n.aM))) 



l<i<V 



e 



\og 13 j 



kj 



log A 



and O* 



log A* 



1/ Using these notations, GAr((fcj)i<j<m) = || ^/f.^^^ ~ ■ O^^,^^ ||^, where | 

j=0 

notes the usual Euclidean norm in R^. Then, with the {£ x £)-identity matrix 



de- 



* 11^ 



j=0 



j=0 
m 

j=0 J 



With P. = ■ (Ll^ ■ L„J-^ ■ L' 



l<i<^ 



from 



< 



mino<j<m(T;+i - T*) N6 



N 



Ell(^r'(*;.*;«)) 



l<i<£ 



J=0 



Now, using the limit theorem 
since k*^^ - k* ~ N5N{r*^i - t* 



N—^oo 



{4 ^(^i' ^i + l))l<i<£ 

oo, and thus 



Af^oo 



AA(0,r(ri,...,r,)) 



(20) 



where ^jy = Op{ipN) as — ^ cxd is written, if for all p > 0, there exists c > 0, such as 
-P^ICa^I < c ■ ipN^ > 1 — p for all sufficiently large N. 

2/ Now, set (fcj)i<j<m € Kytoiv- Therefore, for and N6]y large enough, there exists 
jo £ {1, • • • 5 and (ji, j2) G {1, . . . , m}^ with ji < j2 such that /cj^, < /c*^ — r/Wiv and 
kjo+i > k*^ +r]WN- Thus, 



Let fl* := be the vector such that 



kji kjQ 

h- . — h ■ 



^11^ k*^, - k* 



31 



^ k ■ . — h- 

j=j^ 30 + 1 I^JQ 



/5*exp (a* log(rj aN)) 



kjg^i — k* 
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Then, 



G^((A:,)i<,<^) >|| ^ - (logfi*),<,<, ir + II (logfi*),<,<, 

with Q = (1^-- - (logfir),^^ J' ■ ((logfi*),<^<, - L., ■ 
In the one hand, with (■) defined in (l3|), 

J2-1 



L„„-e^+^ l|2+2g,(2i) 



"-JO 



Oat) + ^ 



A; 



Using the central limit theorems ((61), for and N 6^- large enough, 



E 



< m 



h* — h ■ 



E 



32-1 



j=ji 



E 



(5"^^+' {ri qn) - {ri cat) ) 



+ 



k- , — h- 



E[(^';f^(r,a^)-/5;;(r,a^)°-)' 



E 



"■JO 



- 7^"" 



i2-i 



with 7^ = maxjj{7ff''} (where (7^5'') is the asymptotic covariance of vector Sp^^k^k') 



< C 



Si) 



< 2C 



T]Wn' 



and eq^\k,k')) and C > not depending on A^. Therefore, for large enough, for all 

2 = 1,...,£, 

E[{log{st-^\naM))-\og{Q*)y 
Then we deduce with Markov Inequality that 

ii>^r-(iogi^*),...„ \\'=o, 

>From the other hand, 



' l<i<i 



(22) 



^Jo+^ II 2_ 



Yl ( ("jo log(^i + log A?o) - log ^* 



i^ogn*) ^^^^^ - La^ - Qj^ 

k* — k- 

Define 71 := ^ ■ [3*^_^, for all p G {0, 1, . . . , - Ji - 1}, Ip 



k* 



l3*^+p-i and 7j2-ii+i 



■ 13*^. Then, using Lemma 1411 one obtains 



2 

inf I ^ (^(alog(ria7v) + log/3) - logl]*j } > C* min lajvl 
1=1 



2(°(2)-"(l)) 
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where C > and a^-^-j = maxj=jj_i^...j2 a*, 0:^2) = ^^^j=ji-i,--;j2, ^ conse- 

quence, for satisfying all possible cases of jo, ji and j2, one obtains 

II (logfi:)^<^<, - La, ■ e'^^' ir> C |a^|2(-in»-:— (23) 

Finally, using Cauchy-Schwarz Inequality, 

Therefore, using fl22l) and (|23|1 . since under assumptions of Theorem 12.11 

f, |2(mm, a*-maxi 
= O Oat ^ ' ' ) , 

rjWN V / 

then 

Q = op( II (logfi*),<,<,-L.,-e:-^^ f ). (24) 
We deduce from relations jH]), (ES), ([23]) and ((241) that 

f( min G';v((fc,)i<,<™) > ||a^|2(--»<-'"^^'<)) -^1. 



Proof of Theorem \2.2 : From Theorem 12.11 it is clear that 

kj kj -p 



□ 



F([kj,E]c[k*,k*]) ^ 1 and ,\ ' ^ 1. (25) 



Now, for j = 0, . . . ,m, (xj)i<j<£ G and < £ < 1, let Aj and -B^ be the events such 
that 



A,:={[^„fc;.]c[A;;,fc;+j}n 



A;' kj 



< e 



and B, := y'^L^i (l^^ - ■ 6*) G n(-oo. 



First, it is obvious that 

P(A,)P(5, I A,) < P(5,) < P(5, I A,) + 1 - P(A,). (26) 
Moreover, from (Jll), 

P(i?, I A,) = F(^{ef \k,,k'j)),^,^, e n(-°0'^«] I 

i=i 

p(Ar(0,r(^)(a*,ri,...,r,)) en(-^' 



— 00, Xi 

N-*oo 

1=1 
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Using (l25l) . it is straightforward that P(Aj) — > 1. Consequently. 



N^oo 



i=l 



and therefore \/^^ ^fyf^— -6*) — ^ A/'(0, r*^-'^(Q;*, ri, . . . , r^)) . Now using again 

y Oat V Af^oo 

( l25l) and Slutsky's Lemma one deduces 



- . e.) ^-^^ M(o, r«.(a,- ... .... ,.)). 



Using the expression of Qj as a hnear apphcation of V- ^ , this achieves the proof of Theorem 
[221 ' □ 
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